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Abstract 

The 0(a 3 ) leading logarithmic contribution to the QED radiator 
in the additive form is considered. The effect of the correction on two- 
fermion physics at LEP1 and LEP2 is evaluated and critically compared 
with the one of next-to-leading 0(a 2 ) corrections. A critical comparison 
with existing estimates for the LEP1 energy range is also performed. 
The 0(a 3 ) leading logarithmic corrections turn out to be relevant in 
view of the experimental precision already reached at LEP1 and foreseen 
at LEP2. 
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QED corrections £1X6 , £IS well known, an essential ingredient of precision 
physics at LEP/SLC. In particular, initial-state photonic radiation (ISR) plays 
a central role in the determination of the energy effectively available in the 
center of mass of the e + e~ reaction. Nowadays, this effect is popularly de- 
scribed by the so called Structure Functions (SF) method, pioneered in |J 
and subsequently developed in |2|, 0, |j, j|, || [7], ||. For processes of the kind 
e + e _ — > '-/Z — > // around the Z° resonance, the effect of ISR is crucial in 
the precise determination of the peak cross section. Going above the Z° peak, 
ISR plays a new role, namely it is responsible for the so called U Z° radiative 
return", i.e. the emission of hard photons such that the energy available in 
the kernel reaction is brought back to the Z° mass. Being the Z° radiative 
return caused by ISR, the most precise knowledge of the hard tail of the IS 
photonic spectrum is of utmost importance. It must be noticed that the Z° 
radiative return enhances the contribution of precisely the hard photons that 
reduce the centre of mass energy to the Z° mass. This means that most of the 
events are characterized by the fact that their invariant mass after ISR is very 
close to the Z° mass, or, in other words, they are "LEPl-like" events. More- 
over, once the resonance is produced there is no more particular enhancement 
of any portion of the spectrum: hence final state radiation (and initial-final 
state interference) gives contributions substantially equivalent to the ones it 
already gives at LEP1. Last, but not least, in the case of s-channel processes 
ISR represents a gauge invariant subset of the full set of QED corrections, and 
therefore represents a meaningful subject of investigation. 

In the case in which only a cut on the invariant mass of the event after ISR 
is considered, the SF method provides a very simple recipe for computing the 
corrected cross section as a one- dimensional integration of the proper kernel 
cross section times the so called "radiator" (or "flux function"), namely as 

rl—Xcut 

a(s) = / dxH(x, s)a ((1 — x)s) , (1) 
Jo 

where the radiator is defined as 

H(x,s)= C ±D(z,8)d(—, 8 ) } (2) 

Jl-x Z \ Z J 

D(x,s) being the electron SF. Equation (1), with the definition (2), is a very 
useful tool when considering semianalytical calculations || devoted to data 
analysis; actually, the availability of simple and accurate analytical formulae 
is mandatory for the development of fast fitting programs. 

The radiators available in the literature can be classified as belonging to 
two groups, according to the kind of solution for the electron SF adopted: ad- 
ditive [Q, |], [| and factorized [|], ||, |7| radiators, respectively. As a matter of 
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fact, the additive radiators are extensively implemented in standard compu- 
tational tools used for the two-fermion data analysis at LEP, for instance ||. 
From now on, we will mainly consider this kind of radiators. 

A typical additive radiator consists of an exponentiated part, taking into 
account soft multi-photon emission, plus finite-order leading-logarithmic (LL) 
corrections accounting for hard collinear bremsstrahlung up to O (a 2 ), namely Jl|, 



H NT (x, s) = A 2 Px^ 1 + hi(x, s) + h 2 (x, s), 
hi(x,s) = --p(2-x), 

h 2 (x, s) = \(3 2 



(2 - x) (31n(l -x) -41nx) -4^ '- + x - 6 



(3 = 2{-)[L-l], 



7T, 

2\ 
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L = \n(s/m 

A 2 = l+g)5 1+ g)% 2 . (3) 

The corrections d>i 2, which within the SF formalism are determined at the 
LL level, can be adjusted to take into account process dependent soft plus 
virtual next-to-leading (NL) contributions.^ For processes of the kind e + e~ — > 
7, Z° — > //, they are known in such a way that the radiator of eq. (3) 
reproduces the exact 0(a 2 ) soft plus virtual perturbative results [p| . 

For such processes, a more accurate form of the additive radiator has also 
been derived [12|, taking into account 0(a 2 ) NL hard-photon corrections, in 



such a way that the full 0(a 2 ) perturbative calculation is reproduced, 
namely 

H B (x, s) = H NT (x, s) - h 2 (x, s) + 8% (1 -x,s), (4) 

where 5 2 is defined in eqs. (3.19) and (3.20) of ref. |I2] . The radiator (^) differs 
from (0) by terms which turn out to be numerically negligible at the Z° peak, 
but can in principle be important far from the resonance because of the well 



known relevance of hard photon radiation. Actually, in |T2[ also some ansatz 
for factorized radiators can be found, which reproduce (§) up to the 0(a 2 ) NL 
corrections and take into account additional higher order contributions (more 
on this later). 

1 Actually, also additional pair radiation gives contributions at the 0(a 2 ); it can be 
accounted for by properly redefining the constant 82 for the soft plus virtual part, and 
adding a proper contribution to the hard part of the radiator (jTol . 
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It has to be noticed that the corrections present in eq. (|) and neglected 
in ([3D are dominated by terms of 0(a 2 L), L being the collinear logarithm. 
Being second-order NL corrections, they are in principle of the same order 
of magnitude as the third-order LL ones. These last corrections are already 
known in the literature at the SF level ||, § and for the factorized radiator Q .0 
By using the additive solution of the SF of 0, it is possible to derive the C(/3 3 ) 
corrections to the QED additive radiator. By combining these results with the 
ones already available, we propose the following form for the QED additive 
radiator, including both NL 0(a 2 ) and 0(/3 3 ) contributions: 

H(x, s) = Asfix 13 ' 1 + h x {x, s) + 8" (1 - x, s) + h 3 (x, s), 



h 3 (x, s) 




27 15 1 



+ —x + 4(1 - -x) (tt 2 - 6 In 2 x + 3Li 2 (x)) 



2 4 2 

(6 3 \ f 17 
7 x) +ln 2 (l - x) (-7H 1 — ar 
x 2 J \ x 2 

—6 In a; (6 — x) + 6 ma; ln(l — x) 3a; 



x 



A 3 = A 2 + (^) 3 ^3, 

where ip( n \z) is the n-th order polygamma function, ip( n \z) = d n ip(z)/dz n , 
i/j(z) = T'(z)/T(z). It is worth noting that computing the radiator at C(/3 3 ) 
requires a redefinition of the normalization in front of the exponentiated term, 
A3, which picks up an (D((3 3 ) contribution, 83, originating from the Gribov- 
Lipatov form factor. 

In the following, a sample of numerical results will be shown and com- 
mented. Only results concerning \i cross sections will be considered, since 
the forward-backward asymmetry requires an analysis beyond the aim of the 
present paper. 

In Fig. 1 the separate effects of the NL 0(a 2 ) and 0(/3 3 ) corrections to the 
radiator are shown for the QED corrected cross section as a function of the s' 
cut defined as s'/s > x cu t- The relative deviations with respect to the cross 
section computed by means of the radiator of eq. (§) are shown for several 
centre of mass energies. Both the 0(f3 3 ) (Fig. la) and NL 0(a 2 ) (Fig. lb) 
corrections amount to a contribution of several 0.1% when the Z° radiative 
return is included, but they tend to compensate one another. When the Z° 

2 The 0((3 4 ) corrections to the SF are known to be at 1 x 10~ 4 level, according to || [(]. [uj . 
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radiative return is excluded, or near the Z° resonance, the NL 0(a 2 ) correc- 
tions are confined at the level of 0.01-0.02%, whereas the 0(f3 3 ) ones remain 
at the level of 0.05-0.1%. It should be noted that in the recent analysis shown 
in [0, when studying the radiative corrections to two-fermion production at 
LEP2, only the effect of NL 0(a 2 ) has been taken into account, leading to 
results which, in the light of the present study, are incomplete. 

Figure 2 shows the relative deviations of the cross section computed with 
the full radiator of eq. @ with respect to the one computed by means of the 
one of eq. (|). In Fig. 2a the effect of the newly proposed additive radiator 
close to the resonance is quoted. It is worth noting that at the Z°-peak the NL 
0(a 2 ) plus 0{(3 3 ) corrections introduce a systematic shift of about —0.07%, 
dominated by the 0{(3 3 ) corrections. This effect has not been included in the 
analyses of precision calculations performed in [HJ, but the present experi- 
mental accuracy |17|| requires that it is carefully taken into account. Going 



beyond the Z° peak (Fig. 2b), the effect of the radiator of eq. (|) amounts 
to about —0.1% when the Z° radiative return is excluded, raising to about 
0.25% when it is included. By combining the information of Figs. 1 and 2, one 
concludes that taking into account the NL 0(a 2 ) corrections but neglecting 
the 0{f3 3 ) ones, leads to theoretical predictions for the two-fermion processes 
that are underestimated by about 1% in the inclusive cases. Again, in view of 
the experimental precision foreseen for the inclusive cross sections, this effect 
has to be taken into account in the theoretical predictions. Moreover, it is 
worth noting that in the case of the hadronic cross section, the effect of the 
full additive radiator of eq. (5) as compared to the one of eq. (3), grows up from 
0.25% (see Fig. 2b) to around 0.4% when including the Z° radiative return. 

We notice that, in the case of LEP1 energies, the results of the present 
analysis concerning NL 0(a 2 ) and 0{j3 3 ) corrections confirm the estimates 
obtained in J7|, with the conclusion that, whereas at that time such effects 
were marginal, at present, and in the light of the continuous progress in the 
reduction of the experimental errors, they are significant. As far as the LEP2 
regime is concerned, the present analysis points out the relevance of IS higher 
order hard photon effects, and completes a first attempt to quantify such 
contributions recently performed in [pl|. 



As already pointed out above, in a catalog of the radiators available 



at the time of the LEP1 workshop is given, and several options are offered, 
according to the particular form of the radiators and the contributions taken 
into account. In particular, the options can be grouped as follows: options (A) 
and (B) refer to the radiators reported in the present paper in eqs. (4) and (3) 
respectively, and are in the additive form; options (C), (D) and (E), on the 
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contrary, describe radiators in the factorized form, and are superseded by the 
radiator described in @, The radiator presented here, eq. (5), is intended 
to supersede options (A) and (B), i.e. the additive radiators at present used 
in some computational tools for the LEP analysis. Given this situation, it is 
worth at this point to examine carefully in which relation the best additive 
and factorized radiator are. To this aim, in Fig. 3 a detailed comparison of 
the radiators in themselves as functions of their argument (Fig. 3a) and of the 
corresponding cross sections (Fig. 3b) as functions of the invariant mass cut 
after ISR is shown. The comparison of the radiators shows that their relative 
difference is within 3 x 10~ 4 for most of the spectrum, raising up to 3 x 10~ 3 
in the very hard tail. Analogously, the comparison of the cross sections shows 
that their relative difference is contained within 2.5 x 10~ 4 for centre of mass 
energies ranging from LEP1 to LEP2, and for 0.01 < x cut < 0.99, i.e. for every 
realistic situation. From the present analysis, one can conclude that the two 
radiators are equivalent from the phenomenological point of view. 

The numerical results shown in the present paper point out that the ad- 
ditive radiators implemented in standard computational tools used for the 
two-fermion data analysis at LEP have to be upgraded both for LEP1 physics, 
where they do not fulfil the present precision requirements, and for LEP2 
physics, where they do not describe appropriately the Z° radiative return. 
This is a first result of the present analysis. Moreover, given the results shown 
in Fig. 3, one can conclude that the upgrade can equivalently be performed 
either by definitely substituting the additive radiators by means of the factor- 
ized one of ref. [f|, 0], or by simply modifying the ones already implemented 
with the addition of a few FORTRAN lines, according to the newly proposed 
radiator of eq. (H). This is a second result of the present analysis. At this 
point, the choice is a matter of taste/convenience.Q 

Summarizing, the 0{(3 3 ) contribution to the QED additive radiator has 
been considered, and a new additive radiator, accounting for NL 0(a 2 ) and 
0(/3 3 ) corrections, is proposed. The effect of the 0{f3 3 ) correction on two- 
fermion cross sections at LEP1 and LEP2 has been investigated and critically 
compared with the effect of NL 0(a 2 ) corrections. The 0{(3 3 ) corrections 
turn out to be relevant in view of the experimental precision already reached 
at LEP1; moreover, together with the NL 0(a 2 ) terms, they represent a non- 
negligible contribution to ISR because of the Z° radiative return and in view of 

3 In the approach of ref. @], the "soft photon singularities" are regularized by means of 
the variance-reduction technique known as "control variates" ; for its actual implementation, 
besides the radiator/structure function one needs also their primitives, possibly in analytical 
form in order to save CPU time; hence, for practical reasons, the implementation of an 
additive radiator is preferable. 



5 



the foreseen experimental precision at LEP2. The newly proposed additive ra- 
diator improves considerably the ones implemented in standard computational 
tools for two-fermion processes at LEP. 
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Figure 1: The relative deviation between the cross sections computed by means 
of the radiators of eq. (5) with the substitution &f (1 — x, s) — > h 2 (x, s) (er 3 in 
a) and of eq. (4) (cr 2n i in b), and the radiator of eq. (3) (<7 2 ), as a function of 
the invariant mass cut. 
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Figure 2: The relative deviation between the cross sections computed by means 
of the radiators of eq. (5) (o/) and the radiator of eq. (3) (cr 2 ), as a function 
of the invariant mass cut. The LEP1 (a) and LEP2 (b) cases. 
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Figure 3: The relative deviation between the radiators of eq. (5) H^ d and of 
[4] Hj FS (Fig. 3a), and the corresponding cross sections (Fig. 3b). 
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